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ABSTRACT 


Equations  representing  realistic  electromagnetics  problems  seldom  yield 
exact  solutions,  and  thus  are  usually  treated  numerically.  In  general,  a 
discretization  procedure  such  as  the  method  of  moments  (also  known  as  the 
weighted  residual  method)  is  used  to  convert  the  original  continuous  equation  to 
a  finite-dimensional  matrix  equation.  A  theory  is  presented  that  demonstrates 
the  relation  between  the  eigenvalue  spectrum  of  the  original,  continuous  opera¬ 
tor  and  the  eigenvalues  of  the  method-of-moments  matrix.  In  addition,  an 
equivalence  between  the  finite  difference  method  and  the  method  of  moments  is 
developed  that  permits  the  theory  to  be  applied  to  finite-difference  equations. 
Examples  involving  differential  and  integral  equations  are  used  to  confirm  the 
theory  and  to  illustrate  the  typical  eigenvalue  spectrum  arising  in  electro¬ 
magnetic  field  problems. 
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1 .  INTRODUCTION 


The  treatment  of  realistic  electromagnetic  scattering  problems  is  usually 
beyond  the  scope  of  exact  analytical  techniques.  Instead,  discretization  proce¬ 
dures  such  as  the  finite  difference  method  (FDM)  ^  ,  the  finite  element  method 
(FEM)^,  and  the  method  of  moments  (MoM) ^  (also  known  as  the  weighted  resi¬ 
dual  method)  are  used  to  convert  the  original  continuous  equation  to  a  finite¬ 
dimensional  matrix  equation.  The  solution  of  the  matrix  equation  may  be 
obtained  from  a  straightforward  application  of  direct  algorithms  3uch  as 
Gaussian  elimination,  unless  the  system  is  of  large  order  or  happens  to  be 

poorly  conditioned.  In  the  latter  situations,  preconditioned  iterative  methods 

[4  1 

may  offer  advantages  over  direct  algorithms  . 

Iterative  algorithms  usually  converge  at  a  rate  that  is  highly  dependent 

[4  5  1 

upon  the  eigenvalue  spectrum  of  the  system  matrix  ’  .  This  eigenvalue  spectrum 

is  related  in  turn  to  the  eigenvalue  spectrum  of  the  original  continuous  opera¬ 
tor,  and  to  the  specific  discretization  procedure  (i.e.,  basis  and  testing  func¬ 
tions)  used  to  construct  the  matrix  equation.  Information  about  the  eigenvalue 
spectrum  enables  the  selection  of  a  convergent  iterative  algorithm,  and  may  be 
helpful  for  identifying  optimal  parameters  to  accelerate  the  convergence  of 
iterative  algorithms. 

This  report  presents  a  theory  delineating  the  relationship  of  the  eigen¬ 
values  of  the  continuous  operator  to  those  of  the  matrix  operator.  If  the 
spectrum  of  the  continuous  operator  is  known,  the  theory  can  be  used  to  judge 
the  effects  of  specific  basis  and  testing  functions  on  the  eigenvalue  spectrum 
and  condition  of  the  matrix.  If  the  spectrum  is  unknown,  the  eigenvalues  of  the 


matrix  can  be  used  to  estimate  the  eigenvalue  spectrum  of  the  continuous  opera¬ 
tor,  which  may  promote  a  better  understanding  of  the  continuous  equations  repre¬ 
senting  a  physical  problem.  For  example,  it  will  be  shown  in  Chapter  4  that  the 
eigenvalue  spectrum  of  the  magnetic-field  integral  equation  (MF1E)  is  fundamen¬ 
tally  different  from  that  of  the  electric-field  integral  equation  (EFIE) ,  which 
helps  to  explain  why  certain  iterative  algorithms  are  successful  at  treating  the 
MFIE  but  not  the  EFIE. 

Chapter  2  presents  the  theory  linking  the  eigenvalues  of  the  continuous 
operator  with  those  of  the  matrix.  This  theory  involves  the  explicit  use  of  the 
MoM  discretization  procedure,  and  may  not  appear  to  be  directly  applicable  to 
matrix  equations  constructed  using  the  FEM  or  FDM  procedures.  Chapter  3 
addresses  the  relationship  between  the  different  discretization  procedures,  and 
shows  that  both  the  FEM  and  the  FDM  can  be  interpreted  in  the  context  of  the 
MoM.  Therefore,  the  theory  developed  for  MoM  matrices  can  also  be  applied  to 
finite-difference  matrices  and  finite-element  matrices.  Chapter  4  presents  a 
variety  of  continuous  operators  whose  eigenvalue  spectrums  are  available  from 
analytical  manipulation.  All  of  the  examples  involve  scalar  equations  with 
discrete  spectrums.  Chapter  5  compares  the  eigenvalues  from  the  examples  of 
Chapter  4  to  numerical  results,  in  order  to  confirm  the  theory  and  illustrate 
some  of  its  applications.  Additional  applications  are  discussed  in  Chapter  6. 


3 


2.  EIGENVALUE  PROJECTION  THEORY 


Consider  a  continuous  linear  operator  equation  of  the  form 

L  f  *  g  (2.1) 

where  f  is  an  unknown  function  to  be  determined  and  g  is  a  specified  excitation. 
In  practice,  the  analytical  inversion  of  Equation  (2.1)  may  not  be  possible. 
However,  an  approximate  solution  may  be  obtained  in  terms  of  the  finite  expan¬ 
sion 

N 

f  »  l  a  B  (2.2) 

*•.  n  n 
n“l 


where  the  basis  functions  B  are  specified  and  the  coefficients  a  are  unknowns 

n  n 

to  be  determined.  The  equation  can  be  enforced  approximately  in  terms  of  the 
testing  functions  T^,  producing  a  discrete  system  of  the  form 
N 

l  an  <  Tm,  L  Bn  >  "  <  V  g>  n  *  2 . N  (2.3) 

n»l 

[3  ] 

where  the  brackets  denote  a  suitably  defined  inner  product  .  This  matrix 
equation  can  be  written 

L  a  -  h  (2.4) 

where  the  elements  of  the  matrix  L  are 


mn 


<  T  ,  L  B  > 


m 


and  the  elements  of  the  right-hand  side  are 


h 

m 


<  T 


m 


(2.5) 


(2.6) 
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The  process  of  expanding  the  unknown  function  f  and  enforcing  Equation  (2.1) 
approximately  is  known  as  the  method  of  moments  (MoM).  The  theory  presented 
here  is  similar  to  that  developed  in  Chapter  7  of  [3]. 

The  eigenvalues  of  the  continuous  operator  L  can  be  found  from  the  solution 
of  the  eigenvalue  equation 

L  e  =■  X  e  (2.7) 

where  e  denotes  an  eigenfunction  and  X  the  corresponding  eigenvalue.  If  the 
above  MoM  discretization  procedure  is  applied  to  the  eigenvalue  equation,  the 
unknown  eigenfunctions  become 

N 

c  a  l  e  B  (2.8) 

n  n 
n*l 

and  Equation  (2.7)  reduces  to 

1  en<T»'LBn>-A  I  ’  1  •  2 . N  (2'9> 

n-1  n-1 

This  discrete  system  can  be  written  as  a  matrix  equation 

L  e  -  X  S  e  (2.10) 


where  the  elements  of  the  matrix  L  are 


L 

mn 


<  T 


m 


B  > 
n 


(2.11) 


and  the  elements  of  the  matrix  S  are 


I 

I 

I 


5 


1 


i 


1 


i 

i 


| 


Assuming  the  matrix  in  nonsingular,  Equation  (2.12)  can  be  rearranged  to  yield 


S_1  L  e  -  X  _e 


(2.13) 


Thus,  the  eigenvalues  X  from  Equation  (2.13)  should  approximate  those  of 
Equation  (2.1).  The  latter  is  a  finite-dimensional  matrix  equation  and  will 
yield  at  most  N  independent  eigenvalues.  The  continuous  operator  will,  in 
general,  yield  an  infinite  number  of  eigenvalues.  However,  in  some  sense  the 
dominant  eigenvalues  of  the  original  operator  should  be  approximated  by  the 
eigenvalues  of  the  matrix  operator  S>  ^  L.  In  addition,  the  corresponding  eigen¬ 
functions  should  be  approximated  by  the  expansion 

N 


e  »  l  e  B 
n  n 
n-1 


(2.14) 


.-1 


where  en  are  the  elements  of  the  eigenvectors  of  the  matrix  _L. 

The  above  theory  is  presented  in  the  context  of  the  MoM  discretization,  and 
may  not  seem  directly  applicable  to  finite-difference  or  finite-element  matrices. 
To  show  that  these  relationships  can  be  applied  to  finite-difference  and  finite- 
element  systems,  Chapter  3  discusses  an  equivalence  between  the  three  discreti¬ 
zation  procedures.  Subsequent  chapters  address  the  verification  of  the  theory 
using  examples  from  electromagnetics. 


I 
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3.  DISCRETIZATION  PROCEDURES  FOR  CONTINUOUS  EQUATIONS 

The  previous  chapter  introduced  the  method  of  moments  (MoM)  discretization 
f  3  ] 

procedure  .  The  MoM  requires  the  explicit  introduction  of  basis  and  testing 
functions  in  order  to  discretize  the  domain  and  range  spaces  of  the  continuous 
operator.  Other  widely-used  discretization  procedures  include  the  finite- 
element  method  (FEM)^  and  finite-difference  method  (FDM)  ^  .  This  chapter 
shows  an  equivalence  between  the  FEM,  FDM,  and  MoM  procedures  for  several  simple 
examples. 

The  FEM  has  its  origin  in  the  works  of  Rayleigh  and  Ritz  (1908)  who  devel¬ 
oped  an  approximate  procedure  for  estimating  the  minimum  value  of  a  quadratic 
functional  and  the  function  that  attained  this  value^.  The  Russian  scientists 
I.  G.  Bubnov  (1913)  and  B.  G,  Galerkin  (1915)  observed  that  the  equations 
obtained  by  the  Ritz  procedure  could  also  be  obtained  directly,  i.e.,  without 
regard  to  the  minimization  of  a  functional^.  The  general  procedure  they  devel 
oped  for  the  approximate  solution  of  any  continuous  operator  equation  is  essen¬ 
tially  the  MoM. 

Because  of  the  relationship  between  the  minimization  of  a  functional  by  the 
Ritz  procedure  and  Galerkin’s  method,  the  FEM  and  the  MoM  are  often  equivalent 
procedures.  The  FDM  is  not  based  on  the  same  concepts,  and  cannot  be  made 
exactly  equivalent  to  the  MoM.  However,  the  FDM  can  be  interpreted  as  an 
approximation  to  the  MoM,  at  least  to  dominant  order.  In  other  words,  basis  and 
testing  functions  can  be  found  that  create  a  MoM  matrix  that  closely  approxi¬ 
mates  the  FDM  matrix.  Specific  examples  are  presented  below  to  demonstrate  this 
equivalence  for  FEM  and  FDM  matrices. 


The  two-dimensional  Helmholtz  wave  equation  can  be  written 


2  2 

+  k2f(x,y)  *  g(x,y)  (3.1) 

3x  3y 

where  g  is  some  given  excitation  and  f  is  the  unknown  function  to  be  determined. 

For  vanishing  Dirichlet  boundary  conditions  (i.e.,  f(x,y)  =  0  on  the  boundary  of 

the  region  of  interest),  the  solution  to  Equation  (3.1)  is  the  function  that 

[2  1 

minimizes  the  functional 

F(f)  "  I  //  {  |Vf  |2  -  k2f 2  +  2gf }  dxdy  (3.2) 

The  basic  idea  behind  the  Ritz  procedure  and  the  FEM  is  to  expand  the  unknown 
function  f  in  terms  of  trial  functions  according  to 

N 

f(x,y)  »  l  aR  Bn  (x,y)  (3.3) 

n=*l 


This  expansion  is  substituted  into  Equation  (3.2),  replacing  the  original  func¬ 
tional  by  one  that  now  depends  on  the  N  unknown  coefficients  a} ,  ,  a^ ,  ...  a^. 


If  the  functional  is  differentiated  with  respect  to  each  of  the  unknown 
coefficients,  and  the  result  equated  to  zero  in  order  to  describe  the  stationary 
point,  the  process  yields  the  discrete  system 

N  2 

l  an  //  t7  Bm  •  7  Bn  -  k  BmBn}  dxdy  =  -  //  Bm  s  dxdy  O-4) 

n*l 

m  =  1 ,  2  ,  . . .  ,  N 


This  can  be  written  as  a  matrix  equation 

L  a  ■  h 


(3.5) 


where  the  elements  of  the  matrix  L  are 
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l  »  f|  (7  B  •  7  B  -  k2  B  B  1  dxdy  (3.6) 

mnJJlm  n  m  nJ  7 

and  Che  elements  of  the  right-hand  side  are 

hm  “  "  //  Bm  «  dxd?  (3*7) 

Consider  a  direct  application  of  the  MoM  to  Equation  (3.1),  again  for  the 
case  of  vanishing  Dirichlet  boundary  conditions.  If  the  function  f  is  replaced 
by  an  expansion  of  the  form  identical  to  that  appearing  in  Equation  (3.3),  and 
testing  functions  Tm(x,y)  are  used  to  enforce  the  equation  approximately,  the 
result  can  be  written 

I,  <■„  II  (Tm  V2B  +  k2  T  B  }  dxdy  -  //  I  g  dxdy  0.8) 

n»l 

m  3  1,  2 ,  ...,  N 

If  integration  by  parts  is  used  to  reduce  the  order  of  the  derivatives  appearing 
in  Equation  (3.8),  the  result  is  the  discrete  system 

N  2 

l  an  //  tVTm  *  7Bn  “  k  TmV  ^  "  “  //  Tm  g  d*d?  (3’9) 

n*l 

m  *  1 ,  2  ,  . . . ,  N 

This  system  is  obviously  a  generalization  of  that  appearing  in  Equation  (3.4). 

If  the  testing  functions  are  taken  to  be  identical  to  the  basis  functions,  this 
discrete  system  is  identical  with  that  produced  by  the  functional  minimization 
procedure. 

The  above  comparison  suggests  that  the  FEM  discretization  procedure  is 
equivalent  to  the  MoM  procedure.  The  FDM  is  not  exactly  equivalent  to  the  MoM, 
but  an  approximate  equivalence  can  be  shown.  For  instance,  consider  the  one¬ 
dimensional  wave  equation 


9 


+  k2f(x)  *  g(x) 
dx 


The  conventional  central-difference  formulas  lead  to  the  matrix  equation 


L  a  -  h 


0.10) 


tl] 


(3.11) 


where 


f  2  2 

'  _  _  + 


am 


V. 


and 


A 

0 


h.  * 


m  *  n 


m  »  n  ±  1 


other 


0.12) 


0.13) 


The  interval  size  A  is  assumed  constant. 


Now,  consider  the  same  equation  in  the  context  of  the  MoM.  If  the  basis 
functions 


B(x)-t(x;x  -A  x  ,  x  +  A) 
n  n  •  n  n 


(3.14) 


and  testing  functions 


Vx)  "  A  p(!c;  Xm  'I*  xm  + 


(3.15) 


are  used  to  discretize  Equation  (3.9)  according  to  the  MoM  (these  functions  are 
defined  in  Figure  3.1),  the  matrix  elements  of  Equation  (3.11)  are  given  by 


Wi 


11 


and 


mn 


A 

0 


3 

,  2 

+  4 

k 

m  ■  n 

*! 

k2 

m  ■  n  ±1 

other 

(3.16) 


h 

m 


,  xm  +  A/2 
—  /  g(x)dx 

4  xn  -  4/2 


0.17) 


NoCe  that  Equations  (3.16)  and  (3.17)  agree  to  dominant  order  with  Equations 

(3.12)  and  (3.13),  i.e.,  as  A  ♦  0.  Since  neither  of  the  above  expressions  are 
exact  except  in  the  limit  of  vanishing  A,  Equations  (3.16)  and  (3.17)  may  be  as 
legitimate  as  the  central-difference  formulas  presented  in  Equations  (3.12)  and 

(3.13) .  For  small  values  of  A  these  matrices  are  virtually  identical. 


Now,  consider  the  two-dimensional  wave  equation 

2  2 

-L-l  +  i_i  +  k2f(x,y)  -  g(x,y)  (3.18) 

3x  *y* 


The  conventional  central-difference  approach  yields  the  discrete  operator  with 

m  ■  n 


^-i-+  k2 

A2 


mn 


A' 

0 


m,  n  represents  adjacent  grid  points  (3.19) 


otherwise 


h 


m 


?m> 


(3.20) 


If  the  MoM  procedure  is  used  with  the  modified  pyramid  basis  functions  defined 
in  Figure  3.2,  i.e. , 


.  t  • 

•  *1 


*  > 
VI 

(  V;t‘| 


#i/V1 


•  *1 


'Mi 


BR(x,y)  -  MP  (xn,  yn,  A) 


(3.21) 


Tn(x,y) 


h  *„  - 1-  *m  +  f>  ym 

A 


A  y  +  A> 

2’  7m  2; 


(3.22) 


the  elements  of  the  discrete  operator  are  given  by 


4_  7k_ 

.2  12 

A 


i_  + a! 
2  12 
a 


m,  n  represent  adjacent  grid  points  (3.23) 


k2/48 


m,  n  represent  adjacent  corner 
grid  points 

otherwise 


,  x„  +  A/2  y-  +  A/2 

h  "  “2  J  J  g(x.y)  d*dy 


m  A  xm  -  A/2  ym  -  A/2 


(3.24) 


Again,  the  discrete  operator  produced  by  the  MoM  agrees  with  the  operator  pro¬ 


duced  by  the  conventional  FDM,  at  least  in  the  limiting  case  A  ♦  0. 


This  chapter  demonstrates  that  the  FEM  and  FDM  discretizations  can  be 


interpreted  in  the  context  of  the  MoM,  and  thus  can  also  be  studied  within  the 


scope  of  the  eigenvalue  projection  theory  of  Chapter  2.  Subsequent  chapters 


will  compare  theoretical  eigenvalues  to  those  of  the  finite-difference  and  MoM 


matrices  in  order  to  verify  this  theory. 


*v!v 


*  \  \\  ■ 
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Figure  3.2  Top  view  of  Che  modified  pyramid  function  MP(xn,  yn>  A).  The  func¬ 
tion  is  piecewise  linear  In  x  and  y,  and  vanishes  outside  the  domain 


x  -  A  <  x  <  x  +  A.  y  -  A  <  y  <  y  +  A-  Values  at  inflection 
n  n  n  n 


points  are  shown. 


4.  SOME  CONTINUOUS  OPERATORS  AND  THEIR  EIGENVALUES 


4.1  Introduction 

This  chapter  presents  several  examples  for  which  the  eigenvalues  of  con¬ 
tinuous  operators  are  readily  available  analytically.  Initially,  differential 
operators  of  the  Laplace  and  Helmholtz  type  are  considered.  These  are  fre¬ 
quently  used  in  connection  with  FDM  or  FEM  discretizations.  Subsequently,  we 
examine  the  operators  of  the  electric-field  integral  equation  (EFIE),  the 
magnetic-field  integral  equation  (MFIE),  and  the  combined-field  integral 
equation  (CFIE)  for  electromagnetic  scattering  from  conducting  bodies.  These 
equations  have  been  used  extensively  to  represent  external  scattering  problems. 
All  of  the  examples  presented  here  are  two-dimensional  in  nature,  i.e.,  they 
represent  idealized  structures  that  are  infinitely  long  or  have  no  variation  in 
the  third  dimension.  The  equations  involved  are  scalar  equations.  For  each 
example,  we  identify  the  operator  and  the  corresponding  eigenfunctions  and 
eigenvalues  according  to  the  equ'ation 

L  e  ■  A  e  (4.1) 

where  L  denotes  the  continuous  operator,  e  an  eigenfunction,  and  X  the 
corresponding  eigenvalue.  The  continuous  operators  under  consideration  in  this 
chapter  all  possess  a  discrete  spectrum,  i.e.,  an  infinite  number  of  discrete 
eigenvalues.  In  Chapter  5,  these  same  examples  are  investigated  numerically  to 
show  the  connection  between  the  continuous  operator  and  the  matrix  operator. 


4.2  Poisson's  equation  for  a  grounded  rectangular  cylinder 


Consider  a  grounded  conducting  cylinder  of  rectangular  shape,  containing 
some  given  distribution  of  time-invariant  electric  charge  density.  The  cross- 
sectional  geometry  is  shown  in  Figure  4.1,  where  a  and  b  are  the  rectangle 
dimensions.  The  electric  potential  $(x,y)  must  satisfy  the  Poisson  equation 

2  2 

-^—4  +  ^—4  -  g(x,y)  0<x<a,  0<y<b  (4.2) 

3x  3y 


where  g(x,y)  represents  the  given  charge  distribution,  scaled  to  give  proper 
units.  Since  the  rectangular  cylinder  is  grounded,  the  potential  must  vanish 
everywhere  on  the  cylinder.  Thus,  the  eigenfunctions  of  the  differential  opera¬ 
tor 


L  (*) 

POISSON 


2  2 
3  »  3* 

«  l  l  ■  M  ■  * 

2  2 
3x  3y 


(4.3) 


are 


mn 


3ln  (-Z)  8ln  (JSJ) 


n,  m  >  0 


(4.4) 


and  it  is  easily  verified  that  the  eigenvalues  for  this  case  are 


A 


mn 


if)2  ♦  (f )2 


(4.5) 


Note  that  the  operator  is  positive  definite,  that  Is,  all  the  eigenvalues  are 
real-valued  and  positive.  Some  of  the  smaller  eigenvalues  are  presented  in 
Table  4.1. 

The  eigenvalues  depend  on  the  dimensions  of  the  cylinder,  and  as  the 
cylinder  is  made  larger  they  tend  to  decrease  in  value.  As  the  cylinder  is  made 


smaller,  they  move  away  from  the  origin.  For  any  size  cylinder,  there  are  an 
infinite  number  of  eigenvalues  with  values  greater  than  any  finite  number.  In  a 
loose  sense,  we  can  say  that  they  tend  to  cluster  at  infinity. 


4.3  Helmholtz'  equation  for  a  grounded  rectangular  cylinder 

If  the  source  terms  in  the  previous  example  varied  with  time  in  a  sinu¬ 
soidal  manner,  the  potential  within  the  grounded  rectangular  cylinder  would 
satisfy  the  wave  equation 

2  2 

— I  +  ‘2’  +  k2i|i(x,y)  *  h(x,y)  (4.6) 

3x  3y 

where  h(x,y)  represents  the  excitation,  scaled  to  proper  units.  Assume  that  the 
wavenumber  'k'  is  real-valued,  and  proportional  to  the  reciprocal  of  the  wave¬ 
length.  The  potential  <J>  must  vanish  on  the  cylinder  walls.  In  this  case,  the 
differential  operator  can  be  written 

LHELMHOLTZ(*)  “  "  ~2  ~  ~2  “  k  *  (4,7) 

3x  3y 

The  eigenfunctions  of  this  operator  are  the  functions 

e  (x,y)  -  sin  sin  (~L)  (4.8) 

mn  v  d  '  a  ' 

The  eigenvalues  are  given  by 

'  (f )2  *  (f  )2  -  k2  (4.9) 

b  a 

Depending  on  the  size  of  the  cylinder  relative  to  the  wavenumber  k,  the  operator 
may  not  be  positive  definite.  In  other  words,  for  small  cylinders  all  the 
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eigenvalues  are  positive-valued  (and  in  fact  lie  far  from  the  origin);  as  the 
cylinder  size  grows  with  respect  to  the  wavelength  the  eigenvalues  migrate  to 
the  left  until  the  lowest  crosses  the  origin  onto  the  negative  real  axis.  Thus, 
the  operator  may  be  indefinite  for  this  problem.  For  a  rectangular  cylinder 
measuring  one  wavelength  by  two  wavelengths  in  dimension,  the  smaller  eigen¬ 
values  are  listed  in  Table  4.2.  (The  operator  is  indefinite  for  the  cylinder 
dimensions  of  Table  4.2.) 


4.4  TM-wave  scattering  from  circular  cylinders 

Consider  the  scattering  of  a  time  harmonic  wave  from  a  circular  conducting 
cylinder  of  radius  'a.'  The  cylinder  geometry  is  shown  in  Figure  4.2.  If  the 
excitation  Involves  only  the  TM  polarization  (the  electric  field  parallel  to  the 
cylinder),  an  appropriate  equation  is  the  electric-field  integral  equation 
(EFIE)t31 

E*nC  (*)  -L™e  (Jz)  /  Jz(*’)  Hq(2)  OcR)  a  d*’  (4.10) 

«0 

where 

R  -  2a  jsin  )  |  (4.11) 


i  nc 

Ez  denotes  the  source  or  "incident"  field,  and  the  parameters  ’k’  and  'n' 
describe  the  medium  surrounding  the  cylinder. 

In  this  case,  the  eigenfunctions  of  the  integral  operator  are 


en(0) 


Jn<fr 


(4.12) 


and  Che  corresponding  eigenvalues  are 


XTM,EFIE 


nllka 

2 


J(ka)  H 
n  n 


(2) 


(ka) 


(4.13) 


where  and  denote  Che  n-th  order  Bessel  and  Hankel  functions,  respectively. 

Note  that  the  eigenvalues  are  complex-valued,  and  lie  in  the  right-half 
complex  plane,  as  illustrated  in  Table  4.3  for  a  cylinder  having  one  wavelength 
circumference.  As  the  cylinder  circumference  increases  relative  to  the  wave¬ 
length,  the  eigenvalues  move  around  in  the  complex  plane  as  illustrated  in 
Figure  4.3.  As  can  be  seen  from  Figure  4.3,  the  eigenvalues  may  pass  through 

the  origin  at  certain  values  of  'ka'  (at  the  zeros  of  the  Bessel  function  J  ). 

n 

For  these  values  of  'ka,’  the  cylinder  is  a  resonant  cavity  and  the  equation 

[7  3  ] 

admits  homogeneous  (source-free)  solutions  ’  J.  The  integral  equation  does  not 
have  a  unique  solution  at  these  values  of  'ka.' 


Note  that  for  very  small  cylinders,  where  ka  +  0, 

x™  -  flSti  [X  -  J  f  *  0  (4.14) 


X™  ~  (if)2n+1  +  j  iii  .0,  n  *  0  (4.15) 

and  thus  we  can  say  that  the  eigenvalues  cluster  at  the  origin  as  the  cylinder 
radius  tends  to  zero.  For  very  large  cylinders,  the  eigenvalues  behave 
according  to 

A™  ~  n  cos^  (ka  -  ^0  -  j  n  cos  (ka  -  -jr)  sin  (ka  -  ^-) ,  (4.16) 

ka  +  « 


Thus,  the  eigenvalues  follow  a  circular  path  in  the  complex  plane  for  asymp¬ 
totically  large  cylinders. 
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TABLE  4.3 


FIRST  TEN  DISTINCT  EIGENVALUES  OF  THE  EFIE  OPERATOR 
FOR  A  ONE  WAVELENGTH  CIRCUMFERENCE  CYLINDER,  TM  POLARIZATION. 


8 


346.495 

-  J 

39.964 

114.592 

+  J 

203.433 

7.813 

+  j 

112.240 

0.226 

+  j 

67.395 

0.004 

+  i 

48.773 

0.0 

+  j 

38.488 

0.0 

+  i 

31.854 

0.0 

+  i 

27.194 

0.0 

+  j 

23.735 

0.0 

+  j 

21.062 

4.5  TE-wave  scattering  from  circular  cylinders 


In  the  previous  example,  the  polarization  of  the  incident  electric  field 

was  assumed  to  be  TM  (parallel  to  the  cylinder  axis).  If  the  polarization  of 

the  incident  wave  is  TE  (perpendicular  to  the  axis  of  the  cylinder),  the  EFIE 
[3  1 

takes  the  form 


(V  "  ?k  *  *  (grad  dlv  +  k2)  (  •  V2)  (kR)  a  d<fr’ 

^  *0 

(4.17) 


where  R  is  defined  in  Equation  (4.11).  In  this  case,  the  eigenfunctions  are 


en4>  " 


(4.18) 


and  the  corresponding  eigenvalues  are 


xTE,EFIEm  rjllka  j  ,  (ka)  R  (2)’  (ka) 


(4.19) 


where  the  primes  denote  differentiation  with  respect  to  the  arguments  of  the 
Bessel  and  Hankel  functions.  Table  4.4  presents  the  dominant  eigenvalues  for  a 
cylinder  having  one  wavelength  circumference. 

Note  that  the  eigenvalues  are  complex-valued  and  lie  in  the  right-half 
plane,  as  was  the  case  for  the  previous  example  involving  the  TM  polarization. 
Again,  they  move  around  in  the  complex  plane  as  the  cylinder  radius  is  increased 
relative  to  the  wavelength,  as  illustrated  in  Figure  4.4.  At  certain  values  of 
'ka, '  an  eigenvalue  may  vanish,  and  the  EFIE  has  no  unique  solution  for  that 
cylinder  geometry.  However,  the  behavior  of  the  eigenvalues  is  somewhat  dif¬ 
ferent  from  the  previous  example.  Specifically,  as  the  cylinder  radius  tends  to 
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TABLE  4.4 


FIRST  TEN  EIGENVALUES  OF  THE  EFIE  OPERATOR 
FOR  A  ONE  WAVELENGTH  CIRCUMFERENCE  CYLINDER,  TE  POLARIZATION. 


114.592 

+  j 

203.43 

62.562 

"  j 

167.30 

26.157 

-  j 

313.54 

1.870 

“  j 

526.05 

0.055 

-  j 

727.42 

0.001 

-  j 

921.87 

0.0 

*  j 

1113.89 

0.0 

-  j 

1304.73 

0.0 

-  j 

1494.90 

0.0 

-  j 

1684.65 

X  *',v 
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XQ  *  0 


\  *  -j 

n 


n  *  0 


(4.20) 

(4.21) 


As  the  cylinder  size  increases  asymptotically,  eigenvalues  follow  a  circular 
trajectory  in  the  complex  plane,  much  as  those  of  the  TM  EFIE  operator. 

However,  the  eigenvalues  of  this  operator  tend  to  cluster  at  infinity  (again, 
speaking  loosely)  and  are  thus  different  in  character  from  those  of  the  TM  EFIE 
operator. 


4.6  Alternate  formulation  for  TE-wave  scattering 


The  previous  example  used  the  EFIE  to  represent  TE-wave  scattering  from  a 
conducting,  circular  cylinder.  The  same  problem  could  be  treated  using  an 


alternative  formulation,  the  magnetic-f ield  integral  equation  (MFIE).  The  MFIE 

[3] 


for  this  problem  is 


2n 


-H*nc(*)  -  <J+)  -  /»<♦•>  y*'>  H0(2)UR>  a  d*' 


(4.22) 


,i  nc 


where  R  is  defined  in  Equation  (4.11)  and  denotes  the  incident  magnetic 


field. 


The  eigenfunctions  for  this  operator  are 


e  ($)  ■  e^ni* 
n 


(4.23) 


and  the  corresponding  eigenvalues  are 


XTE,MFIE_  jnka  j  (fca)  H  .(ka) 
n  i  n  n 


(4.24) 
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Table  4.5  shows  the  eigenvalues  for  a  cylinder  having  one  wavelength  circum¬ 
ference. 


The  character  of  the  MFIE  eigenvalues  is  similar  to  those  of  the  EFIE  in 
certain  respects.  The  eigenvalues  are  complex-valued  and  are  located  in  the 
right-half  complex  plane.  Figure  4.5  shows  the  movement  of  the  dominant  eigen¬ 
values  as  the  cylinder  size  is  increased.  Again,  the  eigenvalues  may  pass 
through  the  origin  at  certain  values  of  'a,'  indicating  cavity  resonances  and 
uniqueness  problems  for  those  cylinders.  Similar  to  the  EFIE  operators,  the 
eigenvalues  follow  a  circular  trajectory  as  the  cylinder  becomes  asymptotically 
large. 


However,  as  the  cylinder  radius  tends  to  zero. 


»0  -  1  -  3*  (l*)2  *  1 


(4.25) 


and 


Xn  ~  2  +  i  2n!(n-l) !  ^2  ) 


ka>2n  1 
*  2 


(4.26) 


Thus,  the  MFIE  eigenvalues  tend  to  cluster  at  0.5  +  jO  and  1  +  jO  in  the  complex 
plane.  This  behavior  is  in  contrast  with  the  behavior  of  the  eigenvalues  of  the 
two  forms  of  the  EFIE,  which  tend  to  cluster  at  the  origin  and  at  infinity, 
respectively. 


4.7  Combined-field  equation  for  TE-wave  scattering 


The  EFIE  and  MFIE  operators  for  TE-wave  scattering  from  a  conducting,  cir¬ 
cular  cylinder  are  considered  above.  An  alternate  formulation  involves  the 


TABLE  4.5 


EIGENVALUES  OF  MFIE  FROM  EQ.  (4.22) 
FOR  A  ONE  WAVELENGTH  CYLINDER. 


0.93899  -  j  0.52893 


±  1 


0.60100  +  j  0.22475 


t  2 


0.45486  +  j  0.03795 


t  3 


0.48596  +  j  0.00173 


t  4 


0.49520  +  j  0.00004 


±  5 


0.49775  +  j  0.0 


t  6 


0.49875  +  j  0.0 


i  7 


0.49923  +  j  0.0 


Table  4.6  shows  the  eigenvalues  for  a  cylinder  having  one  wavelength  circum¬ 
ference,  for  the  case  when  the  parameter  a  has  the  value  0.2.  Figure  4.6  shows 
the  movement  of  several  of  these  eigenvalues  as  a  function  of  cylinder  size. 
Note  that  these  eigenvalues  do  not  pass  through  the  origin,  as  do  those  of  the 
EFIE  and  MFIE  operators.  The  CFIE  has  unique  solutions  for  all  values  of 
cylinder  size,  and  thus  may  be  preferred  to  the  other  formulations. 


TABLE  4.6 


EIGENVALUES  OF  THE  CFIE  OPERATOR 
FROM  EQ.  (4.29)  FOR  A 
IX  CIRCUMFERENCE  CYLINDER,  a  -  0.2 


305.92 

-  j 

118.72 

193.65 

+  j 

34.28 

142.32 

-  j 

51.27 

146.83 

-  j 

104.69 

149.26 

-  1 

145.47 

150.02 

-  i 

184.37 

150.32 

-  j 

222.78 

150.46 

-  j 

260.95 

150.54 

-  j 

298.98 

150.59 

-  j 

336.93 

Im(X) 


Figure  4.6  Trajectory  of  three  eigenvalues  of  the  TE  CFIE  applied  to  circular 
cylinders  as  a  function  of  cylinder  size  ka. 


4.8  Summary 


A  variety  of  examples  of  electromagnetic  field  problems  have  been  investi¬ 
gated  to  illustrate  the  type  of  eigenvalue  spectrum  arising  in  practice.  The 
theory  of  Chapter  2  can  be  used  to  relate  the  eigenvalue  spectrum  of  the  con¬ 
tinuous  operator  equation  to  that  of  the  corresponding  matrix  equation.  Since 
these  examples  can  all  be  treated  analytically,  they  permit  the  theory  to  be 
tested  in  a  controlled  manner.  Chapter  5  will  examine  these  examples  again, 
from  a  numerical  implementation,  and  test  the  theory  of  Chapter  2  in  each  case. 


» 


I 


i 
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5.  EIGENVALUE  SPECTRUM  OF  MATRIX  OPERATORS 


Chapters  2  and  3  presented  a  theory  relating  the  eigenvalues  of  a  con¬ 
tinuous  operator  to  those  of  the  finite-difference,  finite-element,  or  method- 
of-moments  matrix.  To  verify  the  theory,  this  chapter  compares  numerical 
eigenvalues  from  such  matrices  with  eigenvalues  of  continuous  operators  pre¬ 
viously  tabulated  in  Chapter  4. 


The  two-dimensional  Poisson  equation  describing  the  potential  within  a 
grounded  conducting  cylinder  containing  free  charge  density  is  presented  in 
Equation  (4.2).  The  geometry  is  displayed  in  Figure  4.1.  The  eigenvalues  of 
the  continuous  operator  equation  for  a  cylinder  geometry  of  dimension  a=*2 ,  b-1 
are  presented  in  Table  4.1.  Consider  a  matrix  approximation  to  this  operator 
based  on  a  finite-difference  method  (FDM)  using  conventional  central-differencing 
formulas^1 The  same  matrix  could  be  obtained  from  the  MoM  using  the  basis  and 
testing  functions  of  Equations  (3.21)  and  (3.22).  In  either  case,  the  matrix 
elements  have  the  form 


l 

mn 


m  =*  n 


m,  n  represent  adjacent  grid  points  (5.1) 


otherwise 


where  A  is  the  interval  between  grid  points  (centers  of  the  basis/testing  func¬ 
tions  in  the  context  of  the  MoM).  According  to  the  eigenvalue  projection  theory 
of  Chapter  2,  the  eigenvalues  of  the  continuous  operator  should  be  approximated 
by  the  eigenvalues  of  the  matrix 


where  the  elements  of  the  scaling  matrix  are  given  by 


r 


7 

12 


m  *  n 


S 


mn 


1 
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m,n  represent  adjacent  grid  points 


(5.3) 


m,n  represent  adjacent  corner 

grid  points 


V. 


0 


otherwise 


Table  5.1  presents  a  comparison  of  the  smallest  eigenvalues  of  the  FDM  matrix 
and  those  of  the  continuous  operator  for  several  different  levels  of  discretiza¬ 
tion.  The  numerical  eigenvalues  appear  to  be  converging  to  those  of  the  con¬ 
tinuous  operator  as  the  order  of  the  discretization  is  increased. 


The  Helmholtz  equation  describing  a  rectangular  geometry  is  presented  in 
Equation  (4.6)  and  eigenvalues  of  the  operator  are  tabulated  in  Table  4.2.  The 
elements  of  a  modified  FDM  matrix  are  presented  in  Equation  (3.23).  The  ele¬ 
ments  of  the  scaling  matrix  used  for  this  case  are  those  given  in  Equation 
(5.3).  Table  5.2  compares  the  eigenvalues  of  the  FDM  matrix  with  those  of  the 
continuous  operator  from  Table  4.2.  Again,  the  matrix  eigenvalues  appear  to 
converge  toward  the  continuous  eigenvalues  as  the  order  of  the  matrix  is 
increased. 


The  electric-field  integral  equation  (EFIE)  describing  the  scattering  of  a 
TM  wave  from  a  circular  cylinder  is  presented  in  Equation  (4.10).  A  matrix 
approximation  to  this  operator  can  be  based  on  the  use  of  pulse  basis  functions 
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TABLE  5.1 


SMALLEST  SIX  EIGENVALUES  OF  S_1L  COMPARED  TO 
THOSE  OF  THE  POISSON  OPERATOR  FOR  a  *  2,  b  *  1. 


21  x  21 
matrix 

55  x  55 
matrix 

105  x  105 
matrix 

Poisson  Operator 
from  Table  4 . 1 

+  12.98 

+  12.63 

+  12.50 

12.34 

+  21.78 

+  20.64 

+  20.25 

19.74 

+  36.99 

+  34.28 

+  33.32 

32.08 

+  46.70 

+  44.31 

43.31 

41.95 

+  59.00 

+  53.89 

51.93 

49.35 

+  59.00 

+  53.89 

51.93 

49.35 

«* 

it 


TABLE  5.2 


SMALLEST  SEVEN  EIGENVALUES  OF  £  L  COMPARED  TO 
THOSE  OF  THE  HELMHOLTZ  OPERATOR,  FOR  a  •  2,  b  *  1,  k 


21  x  21 

55  x  55 

105  x  105 

Helmholtz  operator 
Table  4.2 

-  26.50 

-  26.85 

-  26.98 

-  27.141 

-  17.70 

-  18.83 

-  19.23 

-  19.739 

-  2.49 

-  5.20 

-  6.16 

-  7.402 

7.23 

4.83 

+  3.83 

2.467 

+  19.53 

14.41 

12.45 

9.870 

19.53 

14.41 

12.45 

9.870 

40.50 


30.66 


26.98 


22.207 
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and  Dirac  delta  testing  functions,  according  to 

B  (t)  -  p(i;  ln  -  j,  l  +  4) 

n  n  i  n  i 

T  (l)  -  5(1  -  l  ) 

m  m 


(5.4) 

(5.5) 


The  basis  and  testing  functions  have  been  defined  in  Figure  3.1. 

For  the  basis  and  testing  functions  defined  in  Equations  (5.4)  and  (5.5), 
the  scaling  matrix  defined  in  Equation  2.12  is  an  identity  matrix.  As  a  result, 
the  eigenvalues  of  the  matrix  should  be  a  direct  approximation  to  those  of  the 
original  EFIE  operator.  Table  5.3  presents  a  comparison  of  the  dominant  eigen¬ 
values  of  the  matrix  operator  for  a  system  of  order  30.  The  matrix  eigenvalues 
closely  approximate  those  of  the  continuous  operator. 


For  the  EFIE  representing  a  TE  wave  incident  upon  a  circular  cylinder 
(Equation  (4.17)),  and  the  MoM  matrix  based  on  the  basis  and  testing  functions 


B  (l)  -  t(Z;  £  -  A,  t  .  1  +  A) 

n  n  n  n 

Tm(4,)  “  A  lm  ~  V  ^m  +  I) 


the  scaling  matrix  contains  the  entries 

r 


m,n  represent  adjacent  cells  in  the  model 


(5.6) 

(5.7) 


Table  5.4  presents  a  comparison  of  the  eigenvalues  of  the  matrix  operator  S  L, 


FIRST  ELEVEN  DISTINCT  EIGENVALUES  OF  THE  EFIE  OPERATOR  COMPARED  TO  THOSE 
OF  THE  MOMENT  METHOD  MATRIX  OPERATOR  FOR  A  ONE  WAVELENGTH 
CIRCUMFERENCE  CYLINDER,  TM  POLARIZATION. 

EFIE  (Eq.  (4.13))  30  x  30  matrix 
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TABLE  5.4 

FIRST  ELEVEN  EIGENVALUES  OF  THE  EFIE  COMPARED  TO  THOSE  OF 
THE  MOMENT-METHOD  _L- MATRIX  AND  THE  PRODUCT  MATRIX  S_1L  FOR 
A  ONE  WAVELENGTH  CIRCUMFERENCE  CYLINDER,  TE  POLARIZATION. 
THE  ORDER  OF  S  AND  L  IS  30. 


EFIE  (Eq.  (4.19))  L-matrlx  S  *L-matrix 
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Che  matrix  operator  L_,  and  the  continuous  operator  of  Equation  (4.17).  The 
scaled  eigenvalues  are  a  better  approximation  to  those  of  the  continuous  opera¬ 
tor  than  are  those  of  the  matrix  L^. 

Several  examples  have  been  presented  that  compare  the  eigenvalues  of  con¬ 
tinuous  operators  to  those  of  the  associated  MoM  or  FDM  matrices.  These 
examples  are  intended  to  confirm  the  theory  presented  in  Chapters  2  and  3,  and 
in  addition  provide  a  framework  for  the  future  study  of  the  eigenvalue  spectrum 
of  other  operator  equations.  In  cases  where  eigenvalues  are  not  available  from 
analytical  means,  those  of  the  matrix  ^1.  can  be  used  to  estimate  the  spectrum 
of  the  original  continuous  operator.  This  procedure  may  be  important  for  the 
study  of  typical  eigenvalue  spectrums,  as  few  continuous  operators  arising  from 
practical  electromagnetics  applications  yield  analytical  expressions  for  eigen- 


6.  APPLICATIONS  OF  THE  EIGENVALUE  PROJECTION  THEORY  TO  MATRIX 
CONDITIONING  AND  THE  CONVERGENCE  OF  ITERATIVE  ALGORITHMS 


The  eigenvalue  projection  theory  presented  in  previous  chapters  has  a 
direct  bearing  on  several  applications  of  current  interest.  One  issue  involves 
the  basis  and  testing  functions  used  to  discretize  a  continuous  equation  and 
their  effect  on  the  condition  number  of  the  resulting  matrix  equation.  If 
possible,  it  is  desirable  to  select  basis  and  testing  functions  in  order  to 
achieve  a  well-conditioned  matrix  representation.  Knowledge  of  the  eigenvalue 
spectrum  and  its  behavior  for  limiting  cases  may  be  also  useful  for  the  iden¬ 
tification  of  regions  of  numerical  instability.  For  instance,  certain  numerical 
formulations  fail  when  applied  to  electrically  small  geometries.  An  additional 


issue  involves  the  choice  of  solution  algorithms:  if  an  iterative  procedure  is 
employed  to  solve  the  matrix  equation,  knowledge  of  the  typical  eigenvalue 
spectrum  can  provide  information  essential  to  the  choice  of  a  convergent 


algorithm. 

The  relation  between  the  eigenvalues  of  the  original  operator  and  the  con¬ 
dition  number  of  the  L_-matrix  obtained  from  the  MoM  procedure  is  delineated  in 
Chapter  2.  Since  the  eigenvalues  of  the  continuous  operator  are  fixed  by  the 
specific  equation  and  geometry  of  interest,  the  only  control  over  the  condition 
number  of  the  L-matrix  is  provided  through  the  choice  of  basis  and  testing  func¬ 


tions.  Since  the  eigenvalues  of  the  product  matrix  S_  L  approximate  those  of 
the  original  operator,  the  condition  number  of  ^  *_L  should  be  independent  of  the 
basis  and  testing  functions.  This  suggests  that  the  matrices  S^  and  L  will  tend 
to  compensate  for  each  other.  If  S_  is  poorly  conditioned,  I.  may  also  be  poorly 


conditioned  in  order  that  _S  *_L  have  a  spectrum  in  agreement  with  that  of  the 
original  operator.  As  a  general  rule,  basis  and  testing  functions  should  be 
selected  in  order  to  produce  a  well-conditioned  scaling  matrix  £.  For  example, 
the  S-matrix  produced  by  the  use  of  subsectional  pulse  basis  functions  and  Dirac 
delta  testing  functions  is  often  an  identity  matrix,  and  is  perfectly  con¬ 
ditioned.  On  the  other  hand,  the  use  of  entire-domain  basis  functions  can  lead 
to  a  very  ill-conditioned  ^-matrix  (see  p.  410  of  [9]).  Entire-domain  functions 
are  often  associated  with  poorly  conditioned  L-matrices. 

Suppose  the  continuous  operator  is  such  that  the  eigenvalue  spectrum 
excited  by  the  right-hand  side  of  the  equation  involves  eigenvalues  having  very 
large  and  small  magnitudes.  Such  an  equation  will  be  difficult  to  treat  numeri¬ 
cally,  as  the  L^-matrix  will  be  poorly  conditioned  regardless  of  the  condition 
number  of  S^.  In  this  type  of  situation,  knowledge  of  the  eigenvalue  spectrum 
may  motivate  the  search  for  an  alternative  formulation.  An  example  involving  a 
poorly  conditioned  system  can  be  drawn  from  the  equations  discussed  in  Chapter  4. 
The  EFIE  used  for  TE-wave  scattering  from  circular  cylinders  is  presented  in 
Equation  (4.17),  and  the  associated  eigenvalues  appear  in  Equation  (4.19).  Note 
that  as  the  cylinder  radius  decreases  with  respect  to  the  wavelength,  the  eigen¬ 
values  spread  in  the  complex  plane  according  to  Equations  (4.20)  and  (4.21).  It 
has  been  observed  that  the  numerical  solution  to  the  TE  EFIE  is  numerically 
unstable  for  electrically  small  cylinders ^ 1 ° ^ .  An  alternative  formulation,  the 
TE  MFIE  presented  in  Equation  (4.22),  remains  stable  as  the  cylinder  size 
decreases.  Note  that  the  TE  MFIE  eigenvalues  do  not  spread  as  the  cylinder  size 
i 3  decreased,  but  cluster  in  reasonable  proximity  to  each  other. 
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tn  recent  years,  a  considerable  effort  has  been  directed  toward  the  imple¬ 
mentation  of  iterative  solution  algorithms  for  treating  certain  electromagnetic 
scattering  problems ^ 1 1 ^ .  A  variety  of  integral  equation  formulations  yield 
matrices  having  special  structure  in  the  form  of  slightly  perturbed  Toeplitz 
symmetries'12,13^.  Iterative  algorithms  have  been  used  to  exploit  this  struc¬ 
ture  and  permit  the  analysis  of  geometries  that  are  electrically  larger  than 
those  solvable  by  other  techniques.  However,  simple  iterative  algorithms  some¬ 
times  diverge  when  applied  to  the  general  complex-valued  matrices  arising  from 

,  „  ..  [14,15] 

electromagnetics  problems 


Information  concerning  the  typical  matrix  eigenvalue  spectrum  encountered 
in  these  systems  can  be  used  to  select  a  suitable  iterative  algorithm  for  the 
problem  at  hand.  For  example,  it  Is  well  known  that  the  Jacobi  algorithm  requires 
the  eigenvalues  of  the  iteration  matrix  to  lie  within  the  unit  circle  centered 
at  1+jO  in  the  complex  plane^.  The  eigenvalues  of  the  TE  MFIE  discussed  in 
Chapter  4  appear  to  satisfy  this  criterion,  at  least  for  the  case  of  a  circular 
cylinder.  However,  the  EFIE  formulations  presented  in  Chapter  4  do  not  satisfy 
this  condition.  Therefore,  it  Is  expected  that  the  Jacobi  algorithm  may  be 
divergent  for  systems  representing  the  EFIE  formulation.  (Note  that  the  itera¬ 
tion  matrix  arising  within  the  Jacobi  algorithm  is  not  necessarily  the  L-matrix 
discussed  above,  but  generally  involves  a  splitting  and  subsequent  scaling  in  an 
attempt  to  shift  the  eigenvalues  into  the  necessary  interval^.)  Other  itera¬ 
tive  algorithms  are  capable  of  treating  systems  with  eigenvalues  distributed 
throughout  the  right-half  complex  plane,  but  require  some  knowledge  concerning 
eigenvalue  location  in  order  to  assure  convergence^16^.  Algorithms  that  permit 
the  user  to  specify  parameters  to  accelerate  the  convergence  rate  can  use 
apriori  knowledge  of  the  eigenvalue  spectrum  to  optimize  these  parameters. 


7.  CONCLUSIONS 


An  eigenvalue  projection  theory  is  presented  that  relates  the  eigenvalue 
spectrum  of  continuous  operators  to  the  eigenvalues  of  the  matrix  operator 
arising  from  a  MoM  or  FEM  discretization.  An  approximate  equivalence  between 
the  FDM  and  MoM  procedures  is  developed  which  permits  the  theory  to  be  applied 
to  FDM  matrices.  Numerical  examples  are  used  to  confirm  the  theory  for  integral 
and  differential  equations,  and  to  illustrate  the  typical  eigenvalue  spectrum 
arising  from  electromagnetics  equations. 

Knowledge  of  the  eigenvalue  spectrum  aids  in  the  selection  of  robust 
problem  formulations  and  discretization  procedures.  If  an  Iterative  algorithm 
is  used  to  obtain  the  solution,  such  knowledge  can  also  be  used  to  optimize  the 
algorithm's  convergence  behavior  to  the  problem  of  Interest. 


! 


! 
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